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Abstract 

We consider the new agegraphic model of dark energy with a varying gravitational 
constant, G, in a non-flat universe. We obtain the equation of state and the deceleration 
parameters for both interacting and noninteracting new agegraphic dark energy. We 
also present the equation of motion determining the evolution behavior of the dark 
energy density with a time variable gravitational constant. Finally, we generalize our 
study to the case of viscous new agegraphic dark energy in the presence of an interaction 
term between both dark components. 
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1 Introduction 



Many cosmological observations, such as SNe la [I], WMAP [2], SDSS [3], Chandra X-ray 
observatory j4], etc., reveal that our universe is undergoing an accelerating expansion. To 
explain this cosmic positive acceleration, mysterious dark energy has been proposed. There 
are several dark energy models which can be distinguished by, for instance, their equation 
of state (EoS) (w = during the evolution of the universe. Although the simplest way 
to explain this behavior is the consideration of a cosmological constant [5], the known fine- 
tuning problem [6] led to the dark energy paradigm. The dynamical nature of dark energy, 
at least in an effective level, can originate from various fields, such is a canonical scalar field 
(quintessence) [7], a phantom field, that is a scalar field with a negative sign of the kinetic 
term j8], or the combination of quintessence and phantom in a unified model named quintom 

An approach to the problem of DE arises from the holographic principle that states that 
the number of degrees of freedom related directly to entropy scales with the enclosing area 
of the system. It was shown by 'tHooft and Susskind [10] that effective local quantum field 
theories greatly overcount degrees of freedom because the entropy scales extensively for an 
effective quantum field theory in a box of size L with UV cut-off A. Attempting to solve this 
problem, Cohen et al. showed [11] that in quantum field theory, short distance cut-off A is 
related to long distance cut-off L due to the limit set by forming a black hole. In other words 
the total energy of the system with size L should not exceed the mass of the same size black 
hole i.e. L 3 p^ < LMl 2 where pa is the quantum zero-point energy density caused by UV 
cutoff A and Mp denotes Planck mass ( M 2 = 1/8ttG). The largest L is required to saturate 
this inequality. Then its holographic energy density is given by = 3c 2 M 2 /L 2 in which c is 
free dimensionless parameter and coefficient 3 is for convenience. More recently a new dark 
energy model, dubbed agegraphic dark energy has been proposed [12] (see also [13]), which 
takes into account the Heisenberg uncertainty relation of quantum mechanics together with 
the gravitational effect in general relativity. Following the line of quantum fluctuations of 
spacetime, Karolyhazy p3] proposed that the distance t in Minkowski spacetime cannot be 
known to a better accuracy than 5t = , where A is a dimensionless constant of order 

unity. Based on Karolyhazy relation, Maziashvili proposed that the energy density of metric 
fluctuations of Minkowski spacetime is given by [15] 

PA ~ ~ IT' C 1 ) 

where t p is the reduced Planck time, and M p is the Planck mass. The agegraphic models 
of dark energy have been examined and constrained by various astronomical observations 

DEI El US 

Since we know neither the nature of dark energy nor the nature of dark matter, a mi- 
crophysical interaction model is not available either. However, pressureless dark matter in 
interaction with holographic dark energy is more than just another model to describe an 
accelerated expansion of the universe. Understanding dark energy is one of the biggest chal- 
lenges to the particle physics of this century. Studying the interaction between the dark 
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energy and ordinary matter will open a possibility of detecting the dark energy. It should 
be pointed out that evidence was recently provided by the Abell Cluster A586 in support of 
the interaction between dark energy and dark matter [21]. However, despite the fact that 
numerous works have been performed till now, there are no strong observational bounds on 
the strength of this interaction [22]. This weakness to set stringent (observational or the- 
oretical) constraints on the strength of the coupling between dark energy and dark matter 
stems from our unawareness of the nature and origin of dark components of the Universe. 
It is therefore more than obvious that further work is needed to this direction. 

Previously, it has been shown that the interacting new agegraphic model of dark energy 
can cross the phantom divide [19]. The phantom energy, if it exists, can cause some peculiar 
phenomena e.g. violates the strong energy condition, p + 3p > 0. This leads us to consider 
phantom energy as an imperfect fluid, implying that the phantom fluid could contain non- 
zero bulk and shear viscosities [23]. The bulk viscosities are negligible for non-relativistic and 
ultra-relativistic fluids but are important for the intermediate cases. In viscous cosmology, 
shear viscosities arise in relation to space anisotropy while the bulk viscosity accounts for the 
space isotropy [23] . Generally, shear viscosities are ignored (as the CMBR does not indicate 
significant anisotropics) and only bulk viscosities are taken into account for the fluids in 
the cosmological context. Moreover, bulk viscosity related to a grand unified theory phase 
transition may lead to an explanation of the accelerated cosmic expansion [25] . 
Although the holographic dark energy model with varying gravitational constant has been 
studied in [261 [27], however, until now, in all the investigated agegraphic dark energy models 
a constant Newtons constant G has been considered. The role of G-variation will be ex- 
pressed through the pure number G'/G = Aq, which will be extracted from observations. 
In particular, observations of Hulse- Taylor binary pulsar B1913 + 16 lead to the estima- 
tion G/G ~ 2±4 x 10 _12 yr _1 [HI H2], while helio-seismological data provide the bound 
— 1.6 x l§~ l2 yr~ l < G/G < [13]. Similarly, Type la supernova observations [1] give 
the best upper bound of the variation of G as — 10~ 11 yr~ 1 < ^ < at redshifts z ~ 0.5 
, while astereoseismological data from the pulsating white dwarf star G117-B15A lead 

G 



to 



< 4.10 x 10 ll yr 1 |45j . See also [16] for various bounds on G/G from observational 

data, noting that all these measurements are valid at relatively low redshifts, i.e z < 3.5. 
Since the limits in G-variation are given for G/G in units yr _1 , and since G/G = HG'/G, 
we can estimate substituting the value of H in yr^ 1 [26J. 

Besides, there have been many proposals in the literature attempting to theoretically jus- 
tified a varying gravitational constant. For example in Brans-Dicke theory the gravitational 
constant is replaced by a scalar field coupling to gravity through a new parameter, and it 
has been generalized to various forms of scalar-tensor theories [17] , leading to a considerably 
broader range of variable-G theories. In addition, justification of a varying Newton's constant 
has been established with the use of conformal invariance and its induced local transforma- 
tions [H] . Finally, a varying G can arise perturbatively through a semiclassical treatment of 
Hilbert-Einstein action [19], non-perturbatively through quantum gravitational approaches 
within the "Hilbert-Einstein truncation" [5U], or through gravitational holography [5T| [5^]. 

In the light of all mentioned above, it becomes obvious that the investigation on the 
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interacting new agegraphic dark energy models with varying gravitational constant is well 
motivated. In this paper, we would like to generalize, following [28], the new agegraphic vis- 
cous dark energy models to the universe with spacial curvature in the presence of interaction 
between the dark matter and dark energy with a varying gravitational constant, G. 



2 New ADE with varying gravitational constant 



Soon after the original ADE model was introduced [12] , an alternative model dubbed " new 
agegraphic dark energy" was proposed by Wei and Cai [29], while the time scale is chosen 
to be the conformal time r\ instead of the age of the universe, which is defined by dt = adr], 
where t is the cosmic time. It is worth noting that the Karolyhazy relation 5t = fitp^t 1 / 3 
was derived for Minkowski spacetime ds 2 = dt 2 — dx 2 [TU [13]. In case of the FRW universe, 
we have ds 2 = dt 2 — a 2 dx 2 = a 2 (dr] 2 — dx 2 ). Thus, it might be more reasonable to choose 
the time scale to be the conformal time rj since it is the causal time in the Penrose diagram 
of the FRW universe. The new ADE contains some new features different from the original 
ADE and overcome some unsatisfactory points. For instance, the original ADE suffers from 
the difficulty to describe the matter-dominated epoch while the new agegraphic dark energy 
resolved this issue [29J. The energy density of the new ADE can be written 



3 ™ 2 



where the conformal time is given by 



J a J Ha 2 

If we write t] to be a definite integral, there will be an integral constant in addition. Thus, 
we have f] = 1/a. Let us again consider a FRW universe with spatial curvature containing 
the new agegraphic dark energy and pressureless matter. The Friedmann equation can be 
written 

tt2 k 8irG . . , A , 

H 2 + - = — (p m + Pd) , (4) 

a o 

where k is the curvature parameter with k — — 1, 0, 1 corresponding to open, flat, and closed 
universes, respectively. A closed universe with a small positive curvature (Qk — 0.01) is 
compatible with observations [30]. If we introduce, as usual, the fractional energy densities 
such as 

8irGp m 87rGp D 
then the Friedmann equation can be written 

n m + n D = i + n k . (6) 
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The fractional energy density of the new agegraphic dark energy can also be written 

,2 



n 



D 



H 2 T] 2 



(7) 



We consider the FRW universe filled with dark energy and dust (dark matter) which evolves 
according to their conservation laws 

p D + 3Hp D (l+w D ) = 0, (8) 
Pm + 3Hp m = 0, (9) 

where wd = Pd/pd is the equation of state parameter of new ADE. Taking the derivative of 
Eq. §2§ with respect to the cosmic time and using Eq. (j7|) we have 



Pd 



-Hp 



D 



G'\ 



2— + 77)- 

na G J 



(10) 



where the prime stands for the derivative with respect to x — In a. Inserting this equation 
in the conservation law (jSJ), we obtain the equation of state parameter 



w D 



1 



— G' 
Q D + — ■ 



'111 



3na v " 3G 

The equation of motion for Op can be obtained by taking the derivative of Eq. The 
result is 



D 



D 



na 



(12) 



where the dot is the derivative with respect to the time. The next step is to calculate j^. 
Taking the derivative of both side of the Friedman equation fll]) with respect to the cosmic 
time t, and using Eqs. ©-© and fTTTj) . we obtain 



2 J5 = 3(Q D 



-n 3 / 2 



na 



-i) 



Ql 



n 



D 



(13) 



Substituting this relation into Eq. ( )23|) . we obtain the evolution behavior of the new age- 
graphic dark energy 



D 



D 



D, 



3 wn D 

na 



n k - ^(i + n k - n D ) 



(14) 



(15) 



For completeness, we give the deceleration parameter 

a E_ 
q ~~lrfP~ ' IP' 
which combined with the Hubble parameter and the dimensionless density parameters form 
a set of useful parameters for the description of the astrophysical observations. Substituting 
Eq. (USD in Eq. <$Tb§ we get 



1 3 tt% 2 
a = 2-^ + ^ 



na 2 2G 



(16) 
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3 Interacting new ADE with varying gravitational con- 
stant 



Next we consider the case where the pressureless dark matter and the new ADE do not 
conserve separately but interact with each other. Given the unknown nature of both dark 
matter and dark energy there is nothing in principle against their mutual interaction and it 
seems very special that these two major components in the universe are entirely indepen- 
dent. Indeed, this possibility is receiving growing attention in the literature [3TJ [32j [33] and 
appears to be compatible with SNIa and CMB data [34]. The total energy density satisfies 
a conservation law 

p + 3H(p + p) = 0, (17) 

where p = p m +PD andp = pn- However, as stated above, both components- the pressureless 
dark matter and the new ADE- are assumed to interact with each other; thus, one may grow 
at the expense of the other. The conservation equations for them read 

p D + 3Hp D (l + w D ) = -Q, (18) 
p m + 3Hp m = Q, (19) 

where Q stands for the interaction term. Following [35] we shall assume for the latter the 
ansatz Q = Tp^, with T > which means that there is an energy transfer from the dark 
energy to dark matter. This expression for the interaction term was first introduced in the 
study of the suitable coupling between a quintessence scalar field and a pressureless cold 
dark matter field [311 E2]- We also assume Y = 3£» 2 (1 + r)H where r = p m j ' p D and b 2 is a 
coupling constant. Therefore, the interaction term Q can be expressed as 

Q = 3b 2 Hp D (l+r), (20) 

where 



Combining Eqs. ( jlOl) . (1201 and ([2B with Eq. (1181 we obtain the equation of state parameter 

2 ft;- ,o(l + n fc ) G' 

If we take following [26] < G'/G < 0.07 and assuming Q D = 0.73 and w 0.01 for the 
present time and n = 4, b = 0.1, we obtain —0.87 < wo~ < 0.85 which is consistent with 
recent observations. We can also find the equation of motion for by taking the derivative 
of Eq. (J7|) . The result is 

Q' D = Q D ( -2-^-— v 7 ^] • (23) 

Th(X 




(3 



where 



2^ = m D 



na 



G 



(1 + Q k -Q D ) 



(24) 



Substituting this relation into Eq. ( 123]) . we obtain the evolution behavior of the interacting 
new agegraphic dark energy with variable gravitational constant 



Q D — ill) 



:i-n D ) [3 — y/n D 

na 



Q k -3b 2 (l + Q k )-^(l + Q k -Q D ] 



(25) 



The deceleration parameter is now given by 



Q 



2 2 



n 3 / 2 n q rii 

+ — - -b 2 n + n k ) - —(i + n k - n D ) 

na 2 2 v w 2G V y 



(26) 



Again takeing < G'/G < 0.07 and assuming = 0.73 and fi^ ~ 0.01 for the present time 
and n = 4, b = 0.1, we obtain —0.46 < g— < 0.45 which is again compatible with recent 
observational data |53|. 



4 Interacting viscous new ADE with varying G 

In this section we would like to generalize our study to the interacting viscous new agegraphic 
dark energy model. In an isotropic and homogeneous FRW universe, the dissipative effects 
arise due to the presence of bulk viscosity in cosmic fluids. The theory of bulk viscosity 
was initially investigated by Eckart [36] and later on pursued by Landau and Lifshitz [37] . 
Dark energy with bulk viscosity has a peculiar property to cause accelerated expansion of 
phantom type in the late evolution of the universe [2U ESI |3H]. It can also alleviate several 
cosmological puzzles like age problem, coincidence problem and phantom crossing. The 
energy-momentum tensor of the viscous fluid is 

= PdU^Uv + Pd{9hv + Ufj,u u ), (27) 

where is the four-velocity vector and 

Pd = Pd - 3#£, (28) 

is the effective pressure of dark energy and £ is the bulk viscosity coefficient. We require 
£ > to get positive entropy production in conformity with second law of thermodynamics 
[40] . The energy conservation equation for interacting viscous dark energy is now given by 

p D + 3H(p D +p D ) = -Q, (29) 

which can be written 

p D + 3Hp D (l + w D ) = 9H 2 £ - Q, (30) 
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Combining Eqs. fllOp . (120 p and (1211) with Eq. (1301) we obtain the equation of state parameter 



w D = -1 + — - b 2K —^ + — + — i. (31) 
3na sZd 3G pd 

If we assume £ = aH~ x p D , where a is a constant parameter, then we get 

2 ^- ,,(1 + S\) G' 



^ = 3a - 1 + 3^^-" 2 ^ + 3G- (32) 



The equation of motion for viscous ADE is obtained as 
Q' D = ftp 



i - n D ) (3 - — j + n h - 36 2 (i + n h ) - ^(i + n k - n D ) + 9an$3) 

\ na J G 

5 Conclusions 

In this work we have investigated the interacting new agegraphic viscous dark energy scenario 
with a varying gravitational constant. We have obtained the equation of state and the 
deceleration parameters for both interacting and noninteracting new agegraphic dark energy. 
By consideing non-interacting and interacting cases we have extracted the exact differential 
equations that determine the evolution of the dark energy density-parameter, where the 
G-variation appears as a coefficient in additional terms. 
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